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Abstract 

Toric dynamical systems are known as complex balancing mass action systems in the mathe- 
matical chemistry literature, where many of their remarkable properties have been established. 
They include as special cases all deficiency zero systems and all detailed balancing systems. One 
feature is that the steady state locus of a toric dynamical system is a toric variety, which has a 
unique point within each invariant polyhedron. We develop the basic theory of toric dynamical 
systems in the context of computational algebraic geometry and show that the associated moduli 
space is also a toric variety. It is conjectured that the complex balancing state is a global attrac- 
tor. We prove this for detailed balancing systems whose invariant polyhedron is two-dimensional 
and bounded. 

This paper is dedicated to the memory of Karin Gatermann (1961-2005). 
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1. Introduction 



Toric dynamical systems describe mass-action kinetics with complex balancing states. 
These s ystems have been studied e xtensively in mat hem atical chemistry, st arting with the 
work of iHorn and Jackson! (1 19721) .iFeinberd d 19721) an dlHorn (Il972l . Il973l). and continu- 
ing with the deficiency theory in ( Feinbera . Il979l Il987l Il989l Il995 ) . Mass-action kinetics 
has a wide range of applicati ons in the physical sciences, and now it is beginn i ng to 
play a role in systems biology ( Craciun. Tang and Feinbera . 120061 : iGnacadia et all 2007 ; 
Gunawardena , 20031 ; Sontaa . 2001 ). Important special c ases of these dynamical systems 



include recombination equations in p opulation genetics (lAkin . 1979) and quadra tic dy- 
namical systems in computer science ()Rabinovich. Sinclair and WigdersonlTl992[ ). 

Karin Gatermann introduced the connection between mass-action kinetics and com- 



putational algebra. Our work drew i nspiration both from her publica tions (jGatermann . 
120011 iGatermann and Huberl . 120021 ; iGatermann and Wolfruni . 2005 ) and from her un- 
published research notes on toric dynamical systems. We wh oleheartedly agree with her 
view that "the advantages of toric varieties are well-known" (jGatermannl . l200ll , page 5). 

We now review the basic set-up. A chemical reaction network is a finite directed graph 
whose vertices are labeled by monomials and whose edges are labeled by parameters. 
The digraph is denoted G — (V,E), with vertex set V — {1,2,..., n} and edge set 
E C {(i,j) G V x V : i ^ j}. The node i of G represents the ith chemical complex and 
is labeled with the monomial 



Here Y = (ytj) is an n x s-matrix of non-negative integers. The unknowns c\, C2, . . . , c s 
represent the concentrations of the s species in the network, and we regard them as 
functions Ci(t) of time t. The monomial labels are the entries in the row vector 



9(c) 



Each directed edge (i,j) € E is labeled by a positive parameter Ky which represents 
the rate constant in the reaction from the i-th chemical complex to the j-th chemical 
complex. Note that if there is an edge from i to j and an edge from j to i then we have 
two unknowns Kij and Kji. Let A K denote the negative of the Laplacian of the digraph 
G. Hence A K is the n x n-matrix whose off-diagonal entries are the Kij and whose row 
sums are zero. Mass-action kinetics specified by the digraph G is the dynamical system 



(1) 



A toric dynamical system is a dynamical system (1) for which the algebraic equations 
9(c) ■ A K = admit a strictly positive solution c* 6 K>o- Such a solution c* is a steady 
state of the system, i.e., the s coordinates of 9(c*) ■ A K ■ Y vanish. The requirement that 
all n coordinates of 9(c*) ■ A K be zero is stronger. The first to study toric dynamical 
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systems, Horn and Jackson ( 1972h . called these systems complex balancing mass action 
systems and called c* a complex balancing steady state. A system (1) being complex 
balancing (i.e., toric) depends on both the digraph G and the rate constants Kij. 

Example 1. Let s — 2, n — 3 and let G be the complete bidirected graph on three 
nodes labeled by c\ , C1C2 and c\. Here the mass- action kinetics system (1) equals 



|(ci,c 2 ) = (c\ 
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This is a toric dynamical system if and only if the following algebraic identity holds: 

(K2lK31+K32K21+K23K3l)(Kl3K23 + K2lKl3 + Ki2K23) = («12K 3 2+Kl3K32+K3lKl2) 2 . (3) 

The equation (3) appears in (Horn, 19731 Equation (3.12)) where it is derived from the 



ne cessary and sufficient conditions for complex balancing in mass-action kinetics given 
by [Horn (|l972h . Our results in Section 2 provide a refinement of these conditions. 



Let us now replace G by the digraph with four edges (1,3), (2, 1), (2,3), (3, 1). This 
corresponds to setting «i 2 = K32 = in (3). We can check that, for this new G, the system 
(1) is not toric for any positive rate constants. Note that G is not strongly connected. □ 

Among all chemical reaction networks, toric dynamic al systems have v ery remarkable 
properties. Some of thes e properties are expl ained in (jFeinberd . Il979h . starting with 
Proposition 5.3; see also (jGunawardenal 120031 Theorem 6.4). We shall review them in 
detail in Sections 2 and 3. From our point of view, the foremost a mong these remark able 
properties is that the set Z of all steady states is a toric variety ( Gatermann . l200ll §3). 
Each trajectory of (1) is confined to a certain invariant polyhedron, known to chemists as 
the stoichiometric compatibility class, which intersects the toric variety Z in precisely one 
point c* . In order to highlight the parallels between t oric d ynamical systems and toric 
models in algebraic statistics ( Pachter and Sturmfelsl. 20051 §1.2), we shall refer to the 
steady state c* as the Birch voint: see (jSturmfeld . 1996 . Theorem 8.20). In Example 1, 
the steady state variety Z is a line through the origin, and the Birch point equals 

C* = COnSt • (K12K32 + K 13 K 32 + K 31 K 12 , Kl3^23 + K21K13 + K 12 K 23 ). 

Here the constant is determined because c\ + C2 is conserved along trajectories of (2). 

This article is organized as follows. In Section 2 we develop the basic theory of toric 
dynamical systems within the context of computational algebraic geometry. For each 
directed graph G we introduce the moduli space of toric dynamical systems on G. This 
space parametrizes all rate constants k for which (1) is toric. In Example 1 this space 
is the hypersurface (3). Theorem 9 states that this moduli space is itself a toric variety 
in a suitable system of coordinates. These coordinates are the maximal non-zero minors 
of the Laplacian of G, and the ir explici t form as positive polynomials in the Kij is given 
by the Matrix-Tree Theorem ( Stanley! . Il999l . §5.6). Our results in Section 2 furnish a 
two-fold justification for attaching the adjective "toric" to chemical reaction networks 
with complex balancing, namely, both the steady state variety and the moduli space are 
toric. In addition, the subvariety of reaction networks with detailed balancing is toric. 
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In Section 3 we introduce the Global Attractor Conjecture which states that the Birch 
point is a global attractor for any toric dynamical system. More precisely, we conjecture 
that all trajectories beginning at strictly positive vectors c° will converge to the Birch 
point c* in the invariant polyhedron o f c°. The conjecture is currently open, even for 
deficiency zero systems (cf. Theorem 9). lDe Leenheer. Angeli and Sontae ( 2007 ) found a 
proof for a class of "monotone" deficiency zero networks where the monomials c Vi involve 
distinct unknowns. We prove the conjecture in Section 5 for toric dynamical systems with 
detailed balancing that evolve in a bounded polygon in s-dimensional space. The algebraic 
theory of detailed balancing systems is developed in Section 4. 



2. Ideals, Varieties and Chemistry 

This section concerns the connection between chemical reac tion network theory and 



toric geometry. We use the language of ideals and varieties as in (jCox. Little and O'Sheal . 
l2007l)T Our refere nce on toric geometry and its relations with computational algebra is 
( Sturmfefl 1996). With regard to the dynam ical system (1), we use the notation from 



( Feinber J Tl979. §5) and ( Gunawardena . 20031 . §3) which has the virtue of separating the 



roles played by the concentrations Ci , the monomials c Vi , and the rate constants Kij . 
To study the dynamical system (1) algebraically, we work in the polynomial ring 

Q[c,k] = Q[{ Cl ,c 2 ,...,c s } U {Kij : S E}], 

and we introduce various ideals in this polynomial ring. First, there is the steady state 
ideal (9(c) • A K ■ Y) which is generated by the s entries of the row vector on the right 
hand side of (1). Second, we consider the ideal (9(c) • ^4^) which is generated by the n 
entries of the row vector 9(c) • A K . The generators of both ideals are linear in the Kij 
but they are usually non-linear in the c$. Next, we define the complex balancing ideal of 
G to be the following ideal quotient whose generators are usually non- linear in the k^: 

C G := ((9(c) -A K ) : ( Cl c 2 ■ ■ ■ c s )°°) . 

We have thus introduced three ideals in Q[c, k\. They are related by the inclusions 

(9(c) -A K -Y) C (9(c) -A K ) C C G . 

If I is any polynomial ideal then we write V(I) for its complex variety. Likewise, we 
define the positive variety I/>o(/) and the non- negative variety V>o(I). They consist of 
all points in V(I) whose coordinates are real and positive or, respectively, non-negative. 
Our algebraic approach to chemical reaction network theory focuses on the study of these 
varieties. The inclusions of ideals above imply the following inclusions of varieties: 

V(C G ) C V((9(c)-A K )) C V((9(c)-A K -Y)). (4) 

The definition of C G by means of saturation implies that the left hand inclusion becomes 
equality when we restrict to the points with all coordinates non-zero. In particular, 

V >0 (C G ) = V >0 ((9(c)-A K )). (5) 

Recall from (|Sturmfelsl . [T996h that a toric ideal is a prime ideal which is generated by 



binomials. We soon will r eplace C G by a subideal T G which is tor ic. This is possible b y 
Proposition 5.3 (ii,iv) in (jFeinbergL Il979| ) or Theorem 6.4 (3) in (jGunawardenal 120031 ). 
which essentially state that V>o(Cg) is a positive toric variety. But let us first examine 
the case when C G is a toric ideal already. 
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Example 2. Suppose that each chemical complex appears in only one reaction, and 
each reaction is bi-directional. Hence n — 2m is even and, after relabeling, we have 
E = {(1,2), (2,1), (3,4), (4,3), . ..,(n-l,n), (n,n-l)}. We start with the binomial ideal 

(*(C) • A K ) = ( K 12 C V1 - K 2 lC V2 , K 34: C y3 - K i3 C V \ ... , K n -l <n C Vn - X ~ K n , n -lC yn ). 

The complex balancing ideal Cg is a saturation of (^(c) • A K ), and it coincides with 
the toric ideal of the extended Cayley matrix in the proof of Theorem 7. There are 
many programs for computing toric ideals. For instance, the methods in (jSturmfelsl . 
Il996l §12. A) are available in maple under the command ToricIdealBasis. Explicitly, 
the complex balancing ideal Cg is generated by all binomials n u +c v + — k u ~c v ~ where 

m 

U2i-i,2i(V2i-i ~ V2i) = v and w 2 t-i,2t + "2i,2i-i = for i = 1,2, . . . ,m. (6) 

i=i 

Eliminating ci, . . . , c s from Cg, we obtain the ideal of all binomials k u+ — k u ~ where 
u G N £ satisfies (6) with v = 0. This is the moduli ideal Mq t o be featured in Theorems 
7 and 9 below. It is a prime binomial ideal of Lawrence type (ISturmfelsl . [l996l §7). □ 



Let us next assume that G = (V, E) is an arbitrary digraph with n nodes which is 
strongly connected. This means that, for any two nodes i and j, there exists a directed 
path from i to j. In this case the matrix A K has rank n — 1, and all its minors of size 
(ri— 1) X (n— 1) are non-zero. The next result gives a formula for these comaximal minors. 

Consider any directed subgraph T of G whose underlying graph is a tree. This means 
that T has n — \ edges and contains no cycle. We write k t for the product of all edge 
labels of the edges in T. This is a squarefree monomial in Q[k]. Let i be one of the nodes 
of G. The directed tree T is called an i-tree if the node i is its unique sink, i.e., all edges 
are directed towards node i. We introduce the following polynomial of degree n — 1: 



Ki = E * T ( 7 ) 

T an 2-tree 

The following result is a restatement of the Matrix-Tree Theorem ( Stanley . 19991 §5.6). 



Proposition 3. Consider a submatrix of A K obtained by deleting the i th row and any one 
of the columns. The signed determinant of this (n— 1) x (n— l)-matrix equals (—l)' n ^ 1 Ki. 

This minor is independent of the choice of columns because the row sums of A K are 
zero. Combining Proposition 3 with a little linear algebra leads to the following corollary: 

Corollary 4. The complex balancing ideal Cg contains the polynomials KiC Vj — KjC Vi . 

We now form the ideal generated by these (™) polynomials and we again saturate with 
respect to c\c 2 ■■■ c s . The resulting ideal Tq will be called the toric balancing ideal: 

T G := ( (Ki6» - K jC yi : 1 < i < j < n) : ( Cl c 2 • • • c s )°°) . 

It is thus natural to consider Tg as an ideal in the polynomial subring 

Q[c,K] = Q[ci,...,c s ,K u ...,K n ] c Q[c,k]. 

The claim that this is a polynomial ring is the content of the following lemma. 
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Lemma 5. The polynomials K\, . . . , K n G Q[n] are algebraically independent over Q. 

Proof. Lot K[ G Q[«i, K2, • ■ • , n n ] denote the polynomial obtained from Ki by substi- 
tuting the new unknown ^ for all . We need only verify that the K- are algebraically 
independent, because an algebraic relation among the Ki would be satisfied by the K[ 
as well. Our polynomials are 

K[ = (number of i-trccs in G) • J^[«t- 

The n squarefree monomials Ilt^i K * (f° r * = 1 • • • n ) are algebraically independent 
because an algebraic dependence among these monomials would specify a dependence 
among l/«i, I/K2, • • ■ , Hence, K' X ,K' 2 , . . . , K' n are algebraically independent. □ 

We new discuss the toric balancing ideal Tq. 

Proposition 6. The toric balancing ideal Tq is a toric ideal in Q[c,K]. Moreover, the 
ideal Tq is generated by the binomials K u + ■ c^ uY ^ — K u ~ ■ c^ uY ^+ where u is any row 
vector in Z™ whose coordinate sum U\ + u-i + ■ ■ ■ + u n is equal to zero. 

Proof. Let A denote the edge-node incidence matrix of the complete directed graph on 
n nodes. Thus A is the x n-matrix whose rows are ej — ej for 1 < i < j < n. We also 
consider the nx (n + s)-matrix ( — Y l n ) . The binomials KiC Vj — KjC Vi which define the 
ideal Tq correspond to the rows of the x (n+s)-matrix A - ( — Y I„ ) , and the binomial 
X u + -c( uY }- —K u ~ ■c lyUY ' l+ corresponds to the row vector U-A-(—Y I„ ), where U is any 
row vector of length (™) such that u = U ■ A. The binomial K u + ■ c^- - K u - ■ c^ uY ^+ 
is a Q[cf x , . . . , c^,Ki, . . . , i4T„]-linear combination of the binomials iQc % — KjC Vi . This 
shows that Tq is the lattice ideal in Q[c, K] associated with the lattice spanned by the 
rows of A • ( —Y I„ ), i.e., there are no monomial zero-divisors modulo Tq. To see that 
Tq is actually a toric ideal, i.e. Tq is prime, it suffices to note that Z" +s modulo the 
lattice spanned by the rows of A • ( — Y I„ ) is free abelian of rank s + 1. Indeed, the 
latter matrix has rank n — 1, and its (n — 1) x (n — l)-minors span the unit ideal in the 
ring of integers Z, because each (n — 1) x (n — l)-minor of A is either +1 or —1. □ 

The variety of Tq is a toric variety in SpecQ[c, K], but we continue to regard it as a 
subvariety of C s x C E (or of SpecQ[c, k]). In this interpretation we have 

V >0 (T G ) = V >0 (C G ) - V >0 ((*(c)-A K )). (8) 

Thus Tq still correctly describes the steady state locus of the toric dynamical system. 
The equation (8) holds because the matrix A K has rank n — 1 over the rational function 
field Q(k), and the vector (Ki,K 2 , . ■ . , K n ) spans its kernel under left multiplication. 
Finally, the following elimination ideal is called the moduli ideal of the digraph G: 

M G = T G n Q[4 (9) 

Here Q[k] is the polynomial ring in only the edge unknowns Kij. The generators of Mq 
are obtained from the generators of Cq by eliminating the unknown concentrations Cj. 
For instance, if G is the complete bidirected graph on c\, C1C2 and c\ as in Example 1 
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then the moduli ideal Mq is the principal ideal generated by K1K3 — K\. This coincides 
with condition (3) because K\ — ^21^31 + ^32^21 + ^23^31, and similarly for K2, K$. 

Suppose now that G is an arbitrary directed graph, and let I be the number of con- 
nected components of G. If one of the components Gi f ails to be strong ly connected, 
then V r >o(C < G i ) is empty and hence V > q(Cg) is empty, by ( Feinber j . Il979l Remark 5.2). 
In that case we define Tq and Mq to be the ideal generated by 1. If each connected 
component Gi of G is strongly connected then we define the toric steady state ideal as 



TV, 



((T Gl +T G2 + -.. + T Gl ) : (dca-'-c) 00 ). 



The moduli ideal Mq is defined as before in (9). The equality in (8) still holds and this 
positive variety is in fact non-empty. Here is the first main result of this section: 

Theorem 7. The equations (1) specify a toric dynamical system if and only if the positive 
vector of rate constants Kij lies in the toric variety V(Mq). In this case, the set of steady 
states of (1) with all Ci > equals the set of positive points on the toric variety V(Tq). 



Proof. The positive variety V>o(7b) consists of all pairs (c, k) where k is a strictly 
positive vector of rate constants and c is a strictly positive solution of the complex 
balancing equations $>(c)-A K = 0. The elimination in (9) corresponds to the map of toric 
varieties V(Tq) — > V(Mq) given by (c, k) t— > k. This map is a dominant morphism 
(by definition of Mq), so its image is Zariski dense in V(Mq). The restriction to real 
positive points, ISo(7g) — * V > o(Mq), is a homomorphism of abelian groups (K>o)* 
whose image is dense, so it is the monomial map specified by a matrix with maximal row 
rank. It follows that this restri ction is surjecti ve, and this proves our first assertion. The 
second assertion follows from (|Feinberd . fl979l Proposition 5.3). □ 



We now justify calling V(Mq) a toric variety by writing Mq explicitly as a toric ideal 



•2/1 



in Q[K]. As before, G is a directed graph with n nodes labeled by monomials 6 
We assume that each connected component G\,Gi,...,G\ of G is strongly connected, 
for otherwise Mq = (1). Let Yi denote the matrix with s rows whose columns are the 
vectors yj where j runs over the nodes of the component Gi . We define the Cayley matrix 





(Y t 


Y 2 ■ 


• yA 




1 








Cay G (y) = 





1 







V 





• 1 / 



This is an (s + 1) x n-matrix. Here 1 and are appropriate row vectors with all entries 1 
and respectively. The term "Cayley matrix" co mes from geometric combinatorics, a nd 
it refers to the Cayley trick in elimination theory ( Huber. Rambau and Santosl . 12000) . 

Let S denote the linear subspace of R s which is spanned by the reaction vectors yj — yi 
where G E. This space is known in chemistry as the stoichiometric subspace. We 

write a = dim(S') for its dimension. The quantity 5 :— n~a — l is known as the deficiency 
of the chemical reaction network G. For instance, 5 = 3 — 1 — 1 = 1 in Example 1. 
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Remark 8. The rank of the Cayley matrix Cay G (Y) equals a + l. Hence the deficiency 5 
of the reaction network coincides with the dimension of the kernel of the Cayley matrix. 

The following theorem is the second main result in this section. 

Theorem 9. The moduli ideal Mq equals the toric ideal of the Cayley matrix Cay G (F), 
i.e. Mq is the ideal in Q[K] generated by all binomials K u — K" where u, v G N" satisfy 
Cay G (F) • (u — v) = 0. The codimension of this toric ideal equals the deficiency S. 

Proof. Let Id^ denote the s x s identity matrix and consider the extended Cayley matrix 

/-Id. Yi Y 2 ■■■ YA 

1 ••• 

1 ••• 


v ••• 1 i 

The toric ideal of this matrix is precisely the toric balancing ideal T G , where the unknowns 
C\, C2, . . . , c s correspond to the first s columns. Deleting these s columns corresponds to 
forming the elimination ideal Mq as in (9). This shows that Mq is the toric ideal of the 
matrix Cay G (F). The dimension of the affinc toric variety V(Mq) in C" is equal to 
a + l = rank(Cay G (F)), and hence its codimension equals the deficiency 5 = n—a — l. □ 

We conclude that V > o(Mq) is a positive toric variety of codimension S in M™ - The 
moment map of toric geometry establishes a natural bijection between V>o(M G ) and the 
interior of the Cayley polytope, which is the convex hull of the columns of Cay G (Y). 

In summary, given any chemical reaction network whose components are strongly 
connected, we have shown that the positive toric variety of the Cayley polytope equals 
the moduli space V r >o(M G ) of toric dynamical systems on G. The deficiency S is precisely 
the codimension of this moduli space. In particular, if the deficiency is zero then the 
Cayley polytope is a simplex and (1) is toric for all rate constants K ir Moreover, the 
positive steady states of a toric dynamical system form a positive toric variety V r >o(T' G ). 

3. The Global Attractor Conjecture and Some Biological Applications 

We now consider a fixed toric dynamical system or, equivalently, a chemical reaction 
network (1) which admits a complex balancing state. The underlying directed graph 
G = (V, E) has n nodes labeled by monomials c vi , c V2 , c y ", and we specify positive 
rate constants by fixing a point k° in the moduli space V>o(M G ). We also fix a strictly 
positive vector c° G M> which represents the initial concentrations of the s species. 
The equations (1) describe the evolution of these concentrations over time. We seek to 
understand the long-term behavior of the trajectory which starts at c°, that is, c(0) = c°. 

Let T g (k°) denote the toric ideal in R[c] obtained from T G by substituting the specific 
rate constants k®j G R>o for the unknowns Kij. Then V>o(T g (k )) coincides with the set 
of all steady states of the toric dynamical system (1). The following result is well-known: 
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Proposition 10. [Existence and Uniqueness of the Birch Point] The affine subspace 
c° + S ofM. s intersects the positive toric variety V>o(?g(k )) in precisely one point c* . 



For a proof and references i n the chemistry literature see I Horn and Jackso"nl ( 19721); a 



different proof can be found in iFeinbere (|l979l Proposition 5.3) or iGunawardenal (|2003l 



Proposition 6.4). We remark that variants of Proposition 10 are ubiquitous across the 
mathematical sciences, and the result has been rediscovered many t imes. In statistics, this 
result is known as Birch's Theorem; see ( Pachter and Sturmfeld . I2005L Theorem 1.10). 



To stress the link with toric models in algebraic statistics we call c* the Birch point of 
the toric dynamical system (1) with starting point c°. 

The right hand side of (1) is always a vector in the stoichiometric subspace S = — 
Hi : (i, j) G E}. Hence the trajectory starting at c° stays in the affine subspace c° + S. In 
fact, concentrations remain non-negative, so the trajectory stays in P := (c° + S')nR> . 
We call P the invariant polyhedron. Chemists use the term stoichiometric compatibility 
class for P. The relative interior of P in c° + S is denoted by P° := (c° + S) D R>o- 

Proposition 11. The Birch point c* is the unique point in the invariant polyhedron P 
for which the transformed entropy function 

s 

E(c) = £(c ■ log(c,) - Cl ■ log(c*) -c l+ d) (10) 
i=i 

is a strict Lyapunov function of the toric dynamical system (1). This means the following: 

(a) For all c £ P we have E(c) > and equality holds if and only if c — c* , 

(b) we have dE(c)/dt < along any trajectory c(t) in P, and 

(c) equality in (b) holds at a point t of any trajectory c(t) in P° if and only if c(t) = c* . 



This pr oposition was p roved by IHorn and Jackson (|!972l) . A different proof can be 



found in f Fcinbcrg, Il979f k see especially Proposition 5.3 and its corollaries; see also 
( Gunawardenal [20031 Theorem 6.4) and the paragra ph before it. We suggest comparing 



this with the proof of (jPachter and Sturmfelsl . [20051 Theorem 1.10) 



Any trajectory of the toric dynamical system (1) which starts in the relatively open 
polyhedron P° = (c° + 5)nl s >0 will stay in the closed polyhedron P = (c° + S)nl?, ; 
actually, it is not hard to show that P° is an invariant set. The main conjecture below 
states that an y such trajec tory converges to the Birch point. This conjecture was first 
formulated by I Horn (1974). A steady state x in P° is called a global attractor if any 



trajectory that begins in P° converges to x. 

Global Attractor Conjecture. For any toric dynamical system (1) and any starting 
point c , the Birch point c* is a global attractor of the invariant set P° = (c + S) nR> . 

An important subclass of toric dynamical systems consists of the chemical reaction 
networks of deficiency zero. If the deficiency S = n — a — I is zero then the moduli ideal 
Mq is the zero ideal, by Theorem 9, and (1) is toric for all choices of rate constants. As 
remarked in the Introduction, the Global Attractor Conjecture is open even for deficiency 
zero systems. Our last section is devoted to partial results on the conjecture. First, 
however, we discuss biological examples which illustrate the concepts developed so far. 
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Example 12. [Networks with trivial moduli] We expect that our toric approach 
will be useful for parametric analyses of chemical reaction netw o rks in systems biology . 
Anal yses of this ki nd include (jKuepfer. Sauer and Parrilol 120071 ) , (jGnacadja et al. , 2007 ) 
and (|Sontag| . 120011 ). Many of the explicit examples we found in the literature have trivial 
toric moduli in the sense that either Mq is the unit ideal or Mq is the zero ideal. 

If Mq = (1) then (1) is never a toric dynamical system regardless of what the Kij 
are. This happens when components of G are not strongly connected. Examples in- 
clude Michaelis-Menten kinetics and the covalent modification cycle in (jGunawardena . 

I03l §5). If Mq = {0} then the network has deficiency zero and (1) is always a 
toric dynamical system, regardles s of what the Examples include the cycle in 



Kuepfer. Sauer and Parrilol 2007 . Equation (9)), the monotone networks in ( De Leenheer. Angeli and Sontasi 



2007). 



120071 ). and the following network which is taken from (jGnacadia et al 

The ligand-receptor-antagonist-trap network has s = 8 species and n = 8 complexes. 
This network G has four reversible reactions which we write in binomial notation: 

«15 • C5C6 - K51 ' Cl, K26 ' C Q C 7 - K 62 • C 2 , K37 • C 7 C 8 - K73 ' C3, K48 • C$C 5 - K 84 • C 4 . (11) 

Here I = 4 and a = 4, so that 6 = 0. In the algebraic notation of Section 2, the toric ideal 
Tq equals the complex balancing ideal (^(c)-A K ) and is generated by the four binomials 
in (11). Eliminating c±, c% y . . . , c% as prescribed by (9) yields the zero ideal Mq = {0}. □ 

Example 13. [DHFR catalysis] Here are some examples from systems biology which 
show a more complicated dynamical behaviour. We consider the reaction network in 
(Craciun, Tang and Feinberg, 2006, Figure 5); this reaction network has several posi- 
tive equilibria for some values of the reaction rate parameters (see Craciun, Tang and 
Feinberg, 2006, Figure 7). This reaction network allows for inflow and outflow of some 
chemical species; in the language of deficiency theory, we say t hat one of the complexes of 
this reaction network is the zero complex fsee lFeinberg ( 19791 )). i.e., one of the vectors yi 
is zero. Note that the group A of reactions in this network has almost the same structure 
as mechanism 6 in (Craciun, Tang and Feinberg, 2006, Table 1), shown below: 



E + 51 ^ ESI, E + S2^ ES2, ESI + S2 
ES1S2 -^E + P, SI ^ 0, 52 ^ 0, P -> 0. 



£5152^^52 + 51, 



(12) 



Like the more complicated DHFR catalysis network, the network (12) also has several 
positive equilibria for some values of the reaction rate parameters. It is easy to compute 
the deficiency of this simpler mechanism: the number of complexes is n = 12 (including 
the zero complex), the number of linkage classes is I = 4 (including the linkage class that 
contains the inflow and outflow reactions for the substrates 51, 52 and the product P), 
and the dimension of its stoichiometric subspace is a = 6. Therefore the deficiency of 
the network (12) is 5 = 12 — 4 — 6 = 2. This network cannot be toric for any choice of 
the constant rates because it is not weakly reversible. If we make all reactions reversible 
in (12), then the complexes, the linkage classes, and the stoichiometric subspace do not 
change, so the deficiency of the reversible version of (12) is also 2. 



Example 14. [Recombination on the 3-cube] In population genetics (Akin, 19791 
19821 ). the evolution of a population is modeled by a dynamical system whose right 
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hand side is the sum of three terms, corresponding to mutation, selection and recombi- 
nation. The contribution made by recombination alone is a quadratic dynamical system 
( Rabinovich. Sinclair and Wigderson . 1992f) which can be written in the form (1). In our 



view, toric dynamical systems are particularly well-suited to model recombination. Here 
we consider a population of three-locus diploids, so the underlying genotope of the hap- 
loid gametes is the standard 3-dimensional cube ( Beerenwinkel. Pachter and Sturmfelsl . 



20071 Example 3.9). The eight vertices of the cube are the genotypes. They now play the 



role of the species in chemistry: 

genotypes [000] [001] [010] [011] [100] [101] [110] [111] 

s = 8 

frequencies c\ C2 c 3 C4 C5 c§ C7 eg. 

The recombination network G has n — 16 nodes corresponding to the pairs of genotypes 
which are not adjacent on the cube. There are twelve bidirectional edges, representing in- 
teractions, and we label them using the notation of ( Beerenwinkel. Pachter and Sturmfelsl . 



20071 Example 3.9). Six of the interactions correspond to conditional epistasis: 



[000] H 


- [110] <-> 


[010] - 


- [100] 


«1,2 


C1C7 


- K2,l ■ 


C3C5 


A'i 


= «2,1 


and 


= «1,2 


[001] - 


- [Ill] <-> 


[011] - 


- [101] 


K 3 ,4 


C2C8 


~ K 4 , 3 • 


C4C6 


K 3 


= K 4 ,3 


and i^4 


= «3,4 


[000] H 


- [101] <-> 


[001] H 


- [100] 


K 5,6 


CiC 6 


— K 6,5 • 


C2C5 


K 5 


= «6 : 5 


and liTg 


= K 5,6 


[010] - 


- [Ill] 


[011] - 


- [no] 


«7,8 


C3C8 


- KS.7 ■ 


C4C7 


K 7 


= K 8j7 


and K% 


= K 7 , 8 


[000] - 


- [011] 


[001] - 


- [010] 


«9,10 


C1C4 


— ^10,9 


■ C 2 C 3 


K 9 = 


: Kl0,9 


and If 10 


= K 9,10 


[100] - 


- [Ill] <-> 


[101] H 


- [no] 


Kn,i2 


C5C8 


- «12,11 ' C 6 C 7 


K u = 


K12, 11 


and _R"i2 


= «11,12 



Secondly, we have marginal epistasis, giving rise to the six pairwisc interactions among 

four complexes [000] + [111] [001] + [110] [010] + [101] [100] + [011] 
four monomials K\ 3 -c\c% Ku ■ C2C7 K15 ■ c 3 ce Ki$ ■ C4C5. 

Here if 13, -K14, if 15, ifj.6 are cubic polynomials with 16 terms indexed by trees as in (7). 
By Proposition 3, they are the 3x3 minors of the Laplacian of the complete graph K4: 

' Kl3,14 + tl3, 15 + Kl3,16 — ^13,14 — K 13 : 15 — K 13 : 16 ^ 

— «14,13 Kl4,13 + K 14,15 + K 14,16 — «14,15 — ^14,16 

— Kl5,13 — ^15,14 Kl5,13 + K 15,14 + K 15,16 _ «15, 16 

\ — K 16,13 — ^16,14 — K 16 : 15 Ki6,l3 + Kl6,14 + Ki6 : l5y 

The recombination network G has I — 7 connected components and its deficiency is S = 5, 
as there n = 16 complexes, and the stoichiometric subspace 5 has dimension a = 4. The 
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moduli ideal Mq is minimally generated by 18 binomials. Twelve of them are cubics: 

K S K XX K X5 - K 7 K 12 K W K 6 K 9 K 15 - K 5 K W K W K 4 K XX K X4 - K 3 K X2 K 16 

K 2 K 9 K xi - KiK w Ki 6 K 4 K 7 K U - K 3 K S K 15 K 2 K 5 K U - K X K 6 K X5 

K e K 12 K 13 - K 5 K 1X K X4 K 2 K 12 K X3 - K X K XX K X5 K S K W K 13 - K 7 K 9 K Xi 

K 4 K W K X3 - K 3 K 9 K 15 K 2 K S K X3 - K X K 7 K 16 K 4 K^K X3 - K 3 K 5 K W . 
The remaining six generators of Mq are quartics: 

K 9 K XX K U K X5 - K W K 12 K 13 K X6 K 6 K S K X3 K 15 - K 5 K 7 K U K W 
K 2 K±K 13 K U - K X K 3 K X5 K W K 5 K S K X0 K 1X - K 6 K 7 K 9 K 12 

K X K 4 K W K XX - K 2 K 3 K 9 K X2 K X K A K 6 K 7 ~ K 2 K 3 K 5 K 8 . 

The moduli space (of toric dynamical systems on G) is the toric variety V(Mg) defined 
by these 18 binomials. It has codimcnsion 5 and degree 56. For any recombination rates 
k° 6 V>o(Mg) and any starting point c° in the population simplex A 7 , the trajectory of 
the toric dynamical system (1) stays in the 4-di mensional po lytope (c + S) D A 7 and 
is conjectured to converge to the Birch point c* . I Akinl (|l979f ) calls c* the Wright point. 
It generalizes the classical Hardy-Weinberg equilibrium in the 2-locus system. □ 



4. Detailed Balancing Systems 

In this section we discuss an important subclass of toric dynamical systems called 
detailed balancing systems. Here, every edge of the digraph G exists in both directions. 
We can thus identify G — (V, E) with the underlying undirected graph G = (V, E), where 
E = {{i,j} : (i, j) G -E}. For each undirected edge {i, j} £ E of the graph G we define 
an n x n-matrix A^ 1 '^ as follows. In rows i, j and columns i, j the matrix A^ 1 '^ equals 




and all other entries of the matrix A^ ,J ' are 0. The Laplacian of G decomposes as 

a k = J2 A * ,j} - ( 13 ) 

{i,i}es 

A detailed balancing system is a dynamical system (1) for which the algebraic equations 
^(c) • Ak' 3 ^ = for {i, j} e E admit a strictly positive solution c* 6 R> . In light of 
(13), every detailed balancing system is a toric dynamical system, so the positive solution 
c* is unique and coincides with the Birch point. As it is for toric dynamical systems, the 
condition of being detailed balancing depends on the graph G and the constants Kij . 

We rewrite this condition in terms of binomials in Q[c, «]. The two non-zero entries of 
the row vector ^(c) • A^ 1 '^ are KijC Vi — KjiC Vj and its negative. Moreover, we find 

*( C ) ■ a£J> ■ y = ( KijC y* - K iiC «« ) • ( Vj - yi ), 
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and hence the right hand side of the dynamical system (1) can be rewritten as follows: 

*(c)-A B .y = £ *( c )-4W>.y = 2 {KijdH-Kji^-ivj-vi). (u) 

For a detailed balancing system, each summand in (14) vanishes at the Birch point c*. 

Example 15. We revisit Example 1. Let s = 2,n = 3 and G the complete graph on 
three nodes labeled by cf , cic 2 and c|. The dynamical system (2) is now written as 

^(Cl, C 2 ) = (/tl2cf-K2lClC 2 )-(-l, 1) + (Ki3cf-/t3lc|)-(-2, 2) + (k 2 3CiC 2 -K 3 2C 2 )-(-1, 1). 

This is a detailed balancing system if and only if the following algebraic identities hold: 

K? 2 «31 - «21 K 13 = «23 K 31 - «32«13 = «12«32 ~ K21«23 = 0. (15) 

This defines a toric variety of codimension two which lies in the hypersurface (3). □ 

To fit our discussion into the algebraic framework of Section 2, we now propose the 
following definitions. The detailed balancing ideal is the following toric ideal in Q[k, c}: 

Tq := ( (nijC Vi — KjiC Vj | {i,j} G E) : ( Cl c 2 • • • c s )°° ) . (16) 

The corresponding elimination ideal in Q[k] will be called the detailed moduli ideal: 

M G := f G n Q[k]. 

The ideal Tq is toric, by the same reasoning as in Proposition 6. The detailed moduli 
ideal Mq is a toric ideal of Lawrence type, as was the ideal in Example 2. Note, however, 
that the ideals Tq and Mq are toric in the original coordinates Kij. Here, we did not 
need the transformation to the new coordinates K\, . . . , K n in (7). 

Using the ring inclusion Q[K, c] C Q[k,c], we have the following inclusions of ideals: 

T G C T g and M G C Mq. 

Here the equality holds precisely in the situation of Example 2, namely, when each 
chemical complex appears in only one reaction and each reaction is reversible. In general, 
as seen in Example 15, the corresponding inclusion of moduli spaces will be strict: 

V >Q (M G ) c V >0 (M g ). 

In words: every detailed balancing system is a toric dynamical system but not vice versa. 

The following characterization of detailed balancing systems will be used in the next 
section. If L is any vector in R s and c the unknown concentration vector then we write 

L * c := (iici, L 2 c 2 , . . . , L s c s ). 

Lemma 16. A toric dynamical system is detailed balancing if and only if all the binomials 
KijC Vi —KjiC Vj in (16) have the form (L*c) Vi ~ [L*c) Vj , for some positive vector L G K> - 
Thus, a detailed balancing system is a toric dynamical system of the special form 

dr 

- = Yl ((L*c)^-(L*c)^).( yj - yi ). (17) 
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Proof. The if-direction is easy: if our binomials have the special form (L * c) Vi — (L * c) Vj 

then c* = (l/Li, 1/L 2 , ■ . ■ , 1/7 S ) is a positive solution to the equations \I/(c) • A^ 1 '^ = 0. 
Conversely, for the only-if direction, we define L as the reciprocal of the Birch point 
L = (1/c*, 1/c*., . . . , 1/c*), and the result follows the fact that c Vl ~ V] = [c*) Vi ~ y i remains 
valid for all stationary points c of the system (1) as the starting point c(0) varies. □ 

We now fix a detailed balancing system (17) with a particular starting point c(0). 
Then the trajectory c(t) evolves inside the invariant polyhedron P = (c(0) + S) n K> . 
Consider any acyclic orientation E' C E of the graph G. This means that -B' contains 
one from each pair of directed edges and (j, i) in E, in such a way that the resulting 
directed subgraph of G has no directed cycles. The acyclic orientation E' specifies a 
stratum S inside the relatively open polyhedron P° = (c(0) + S) (1 R> as follows: 

S := {c€ P° \ (L* c) m > (L* c) Vi for all in E'}. 

The invariant polyhedron P is partitioned into such strata and their boundaries. We 
are interested in how the strata meet the boundary of P. Each face of P has the form 
Fi := {c G P | Cj = for i G 7} where 7 is subset of {1,2,..., s}. This includes i*0 = P. 

Lemma 17. Consider a detailed balancing system (17) and fix an acyclic orientation E' 
of the graph G. If the closure of the stratum S corresponding to E' intersects the relative 
interior of a face Fi of the invariant polyhedron P, then there exists a strictly positive 
vector a G R> such that ^2 keI (yjk — Vik) • ctk > for all directed edges (i,j) in E' . 

Proof. We proceed by contradiction: assume that the inequalities ^2 keI (yjk —yik)ctk > 
have no strictly positive solution a G K> - By Linear Programming Duality (Farkas' 
Lemma), there is a non-negative linear combination v = j)eE' ^ij(Vj — Vi) sucn * na t 
the following two conditions on v hold: (a) supp(w + ) PI I = 0, and (b) supp(v~) contains 
some j G I. We shall prove the following two claims, which give the desired contradiction: 

Claim One: If c is a point in the relative interior of Fj , then (L * c) v+ > (L * c) v ~ . 

Since (L * c)i — if and only if i G I, and (L * c)j > for all j £ I, (a) implies that 
[L * c) v+ is strictly positive, while (b) implies that (L * c) v ~ — 0, and we are done. 

Claim Two: If c is a point in the closure of the stratum S, then (L*c) v+ < (L*c) v - . 

Consider any point s G S. By the construction of v, the following equation holds: 

(L*s) v = {L * s)^^ E ' Xiiivi ~ Vi) = Y[ ((L*s) v *- yi ) Xii . (18) 

Recall that (L * s) Vj ~ Vi < 1 for each oriented edge G E'. Also, each Ajj is non- 

negative, so ((L* s) y i~ Vi ) Xi i < 1. Using (18), this implies that (L*s) v < 1, and therefore 
(L * s) v+ < (L * s) v - . By continuity we can replace s by c in this last inequality. □ 

The vector a G K> in Lemma 17 will play a special role in the next section. In Corol- 
lary 18 below we regard a as a vector in R> by setting aj = for all j G {1, . . . , s}\I. 

Corollary 18. Let c(t) be a trajectory of a detailed balancing system (17) on the invari- 
ant polyhedron P, and suppose that a point c(to) on this trajectory lies both in the closure 
of a stratum S and in the relative interior of a face Fj of P. Let a G K> be the vector 
obtained as in Lemma 17. Then, the inner product (a, ^f(to)) * s non-negative. 
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Proof. Let E' denote the orientation which specifies S. The velocity vector ^| (to) equals 
]T ((L*c(t ))^-(L*c(t ))^)-( yj - yi ). 

Since c(to) is in the closure of the stratum S, we have (L * c(t )) yi — (L * c(to)) Vj > 0. 
We also have (a, yj — yi) > because a comes from Lemma 17. This implies 

<«, §(*o)> = E ((i* C (io)r"(i*c(M) w )-(« ; »i-Wi> > 0. 

This is the claimed inequality. It will be used in the proof of Theorem 23. □ 



5. Partial Results on the Global Attractor Conjecture 

This section contains what we presently know about the Global Attractor Conjecture 
which was stated in Section 3. This conjecture is proved for detailed balancing systems 
whose invariant polyhedron is bounded and of dimension two. We begin with some general 
facts on trajectories of toric dynamical systems, which are interesting in their own right. 

Consider a fixed toric dynamical system (1) with strictly positive starting point c(0) = 



c 



€ K> - The trajectory c(t) remains in the invariant polyhedron P = (c + S) fl . 



<->o- 



Recall that any face of P has the form Fi := {c G P | Cj = if i £ /}, where / C {1, . . . s}. 
The boundary dP of P is the union of all faces Fj where I is a proper subset of {1, . . . , s}. 
For positive e, the e-neighborhood in P of the boundary of P will be denoted by V £ (dP). 

We note that the transformed entropy function (10) can be extended continuously to 
the boundary of P, because c^logc^ — * as Cj — » + . Equivalent formulations of th e 



following result are well known. For instance, see lSiegel and Chenl (|1994I ); ISontad (|2001r ). 



Proposition 19. Suppose that the invariant polyhedron P is bounded and the distance 
between the boundary of P and the set {c(t) G P \ t > 0} is strictly positive. Then the 
trajectory c(t) converges to the Birch point c* of P. 



Proof. We assume that c(t) does not converge to c*. Let e > be such that c(t) ^ 
V £ (dP) for all t > to. The strict Lyapunov function (10) ensures that there exists a 
neighborhood V £ >(c*) of the Birch point c* such that all trajectories that visit V e '(c*) 
converge to c*. Then c(t) £ V e i (c*) for all t > to- Denote the complement of the two open 
neighborhoods by Po := P \ (V e (dP) U V e > (c*)). Then the non-positive and continuous 
function c i— > (VE- )(c) does not vanish on P by Proposition 11, so it is bounded above 
by some —5 < on Po. Therefore, the value of E(c(t)) decreases at a rate of at least 8 
for all t > to, which implies that E is unbounded on P - This is a contradiction. □ 

Given a trajectory c(t) of (1), a point c £ P is called an uj-limit point if there exists a 
sequence t n — > oo with linin^oo c(t n ) = c. Proposition 19 says that if the trajectory c{t) 
does not have any w-limit points on the boundary of P, then it must converge to the 
Birch point c* . Thus, in order to prove the Global Attractor Conjecture, it would suffice 
to show that no boundary point of P is an cj-limit point. We first rule out the vertices. 
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Proposition 20. Let r be a vertex of P and consider any e > 0. Then, there exists a 
neighborhood W of r such that any trajectory c(t) with starting point c(0) = c° satisfying 
dist(c , r) > e, does not visit W for any t > 0. 



Proof. The following set is the intersection of a closed cone with a sphere of radius one: 

V := I -At I v e S\{0\ and r + v lies in p\ . 
UMI J 

Hence V is compact. We set I — {j £ {l,...,s} : Tj = 0}. For each » £ V, the ray 
7„(t) := r + tv extends from the vertex 7„(0) = r into the polyhedron P for small t > 0. 
We consider how the transformed entropy function changes along such a ray: 

d s 

j t E(iv{t)) = 5>,(io g (o + to 3 )) + J^tyiGgfo+tty) - H 1o s( c >j) 

je/ i=i 

= (ZjElVj)-\0g(t) + w(t), 

where the function w(t) admits a universal upper bound for t close to and v E V. 
For each j 6 / we have Uj > because 7"j = and r + tv £ P for small t > 0. Also, 
since w points into P, there exist j £ I with Vj > 0. Thus, the function £j* 6 jUj has a 
positive minimum over V. It follows that ^E{^ v {t)) tends to — oo for t — ► 0. There exists 
to < £ such that for all w £ V the function t i— » J5(r + to) decreases for < < < <o- 
So, J5(r) > /i := max„ e v _B(r + tov). On the other hand, E is continuous, so there is 
a neighborhood VF of the vertex r (contained in {r + to | i < to, v S V}) such that 
£J(c) > (P(r) + fj,)/2 for all c G W. Since P decreases along trajectories, we conclude 
that no trajectory c(t) that starts at distance > e from the vertex r can enter W. □ 



Remark 21. Chemical reaction networks for which P is bounded are called conservative. 
For conservative networks, th ere exists a posit ive mass assignment for each species that 
is conserved by all reactions ( Feinbere . 19791 ). On the other hand, if € P, then the 



reaction network is not conservative. Proposition 20 ensures that, for a toric dynamical 
system, complete depletion of all the concentrations C\,C2, ...,c s is never possible. 

Lemma 22. Suppose that P is bounded and that the trajectory c(t) has an w-limit point 
on the boundary of P. Then for any e > there exists a positive number t £ > such that 
c(t) belongs to V £ (dP) for all t > t s . In other words, the trajectory c(t) approaches the 
boundary. 



Proof. Suppose that for some e > there exists a sequence t n — > oo such that c(t n ) ^ 
V e (dP) for all n. As P is bounded, the trajectory c(t) has an w-limit point p S P\V e (dP). 
On the other hand, c(t) also has an w-limit point on the boundary of P. Consider a ball 
B<z${p) of radius 28 around p, whose closure lies fully in the relative interior of P. The 
trajectory c(t) enters and exits the neighborhood B$(p) of p infinitely many times, and 
also enters and exits the neighborhood P\E>2s(p) of the boundary infinitely many times. 
The trajectory c(t) travels repeatedly between these two sets which are at distance 6 
from each other. Note that \dc/dt\ is bounded above, and VE ■ dc/dt is bounded away 
from zero on the annulus B2s(p)\B$(p). Then, as in the proof of Proposition 19, each 
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traversal between the neighborhoods decreases the value of E(c(t)) by a positive amount 
that is bounded away from zero. This contradicts the fact that E is bounded on P. □ 

We shall now prove the main result of this section. Admittedly, Theorem 23 has 
three rather restrictive hypotheses, namely, "dimension two," "bounded polyhedron," and 
"detailed balancing." At present we do not know how to remove any of these hypotheses. 

Theorem 23. Consider a detailed balancing system (17) whose stoichiometric subspace 
S = R{yj — yi | (i,j) e E} is two-dimensional and assume that the invariant polygon 
p = (c» + s)n M> is bounded. Then the Birch point c* is a global attractor for P. 

Proof. By Proposition 19, we need only rule out the possibility that the trajectory c(t) 
has an w-limit point on the boundary of P. Proposition 20 gives the existence of open 
neighborhoods of the vertices such that no trajectory c(t) that starts outside them can 
visit them. Let V denote the union of these neighborhoods. Suppose now that c(t) has 
an u-limit point on dP. That limit point lies in the relative interior of some edge F of 
P. Let F e denote the set of points in P which have distance at most e from the edge F. 

We claim that there exists e > and t E > 0, such that the trajectory c(t) remains 
in the subset F £ \V for all t > t £ . This is true because c(t) belongs to the neighborhood 
V e (dP) of the boundary for t 3> 0, by Lemma 22, and hence c(t) belongs to V e (dP)\V for 
t ^> 0. But this implies that c(t) belongs to F E \V for t ^ because F e \V is a connected 
component of V £ (dP)\V for e sufficiently small. This uses the dimension two assumption. 

Consider the closures of all strata S that intersect the relative interior of F. After 
decreasing e if necessary, we may assume that the union of these closures contains the 
set F e \V, which contains the trajectory c(t) for t > t £ . To complete the proof, we will 
show that the distance from c(t) to the edge F never decreases after c(t) enters F e \V. 

Any stratum S whose closure intersects the relative interior of F contributes a vector 
a = a(S) which satisfies the statement of Lemma 17 for F = Fj. The orthogonal 
projection of a(S) into the two-dimensional stoichiometric subspace is a positive multiple 
of the unit inner normal a E S to F in P. By Corollary 18 we have (a(S), > 

and hence (ao, 4f (t)) > for t > t £ . Therefore the distance from c(t) to F cannot 
decrease. This is a contradiction to the assumption that F contains an w-limit point. □ 
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